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Abstract 

Two-magnon Raman scattering is a useful tool to verify recent suggestions 
concerning the value of the interplanar exchange constant in antiferromag- 
netic two-layer systems, such as YBa2Cu30Q+x- We present a theory for 
Raman scattering in a two-layer antiferromagnet. We study the spectra for 
the electronic and magnetic excitations across the charge transfer gap within 
the one-band Hubbard model and derive the matrix elements for the Ra- 
man scattering cross section in a diagrammatic formalism. We analyze the 
effect of the interlayer exchange coupling J2 for the Raman spectra in Aig 
and Big scattering geometries both in the non-resonant regime (when the 
Loudon-Fleury model is valid), and at resonance. We show that within the 
Loudon-Fleury approximation, a nonzero J2 gives rise to a finite signal in 
Aig scattering geometry. Both, in this approximation and at resonance, the 
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intensity in the Aig channel has a peak at small transferred frequency equal 
to twice the gap in the spin-wave spectrum. We compare our results with 
experiments in YBa2Cu30Q,i and Sr2Cu02Cl2 compounds and argue that 
the large value of J2 suggested in a number of recent studies is incompatible 
with Raman experiments in Aig geometry. 
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I. INTRODUCTION 



Since the discovery of high-Tc superconductivity [jT| a lot of work has been done in 
an attempt to understand the pairing mechanism. Most of the existing theories consider 
boson-mediated pairing between electrons in the same Cu02 plane, but there also exist 
arguments in favor of the pairing between electrons in adjacent Cu02 planes These 
arguments are mostly applied to Y Ba2Cu^0Q+x iYBCO) compounds in which the unit 
cell contains pairs of Cu02 planes separated by a charge reservoir. A strong magnetic 
coupling between the planes of a bilayer was also suggested |^-|^ as a possible source for 
an experimentally observed strong downturn renormalization of the low temperature Pauli 
susceptibility of holes in underdoped YBCO compounds ("spin-gap" phenomenon) as well 
as for the maximum in the spin-lattice relaxation rate at temperatures well above 0. 

An essential input parameter for these theories is the value of the interplane hopping 
amplitude or the Cu — Cu superexchange interaction between the two Cu02 planes. In this 
paper, we will argue that the Raman scattering experiments in two-layer compounds allow 
an estimate for the value of the interlayer exchange coupling. 

The Raman scattering in single-layer parent high-Tc compounds has been intensively 
studied over the last few years ||7|-[T^. A large number of studies has been performed to 
understand the Raman spectra in different scattering geometries. In Big geometry (sj = 
(x + y)/\/2, ej = (x — y)/V^ [|l3l, where ejj are polarization unit vectors of the incident 
and scattered photons) the dominant feature of the magnetic Raman intensity profile is 
a peak at about SOOOcm^^, which is attributed to a two-magnon scattering process [jl4|. 
The two-magnon peak has been observed in all parent high-T^ compounds. Besides, the 
experiments also found a strong Raman signal in the Aig geometry (ej = ej = (x + y)/\/2) 
p!5| . The Aig signal has a maximum at about the same frequency as in Big geometry, but 



the width of the peak is larger and its intensity is about a quarter as strong. 

A traditional framework for the understanding of Raman experiments is the Loudon- 
Fleury theory [|T6[ which describes the interaction of light with only spin degrees of freedom. 



3 



This theory explains a peak in Big geometry but predicts that there should be no scattering 
in Aig geometry. Recently, however, it was found that the Loudon-Fleury approach has 
to be modified because Raman experiments are mostly performed near the resonant regime 
where photon frequencies are close to the charge transfer gap of the insulating compounds, 
and one can by no means neglect electronic degrees of freedom. In this regime, the diagrams 
which are neglected in the Loudon-Fleury theory, and which contribute to both Aig and Big 
scattering are actually more important than the diagrams included in the Loudon-Fleury 
theory. 

In this paper, we will study magnetic Raman scattering in two-layer antiferromagnetic 
insulators. We will show that in the presence of the interlayer exchange coupling, J2, the 
Raman scattering profile in the Aig scattering geometry acquires qualitatively new features. 
In particular, a nonzero J2 gives rise to a nonzero Raman intensity in the Aig scattering 
geometry already in the Loudon-Fleury approximation. Moreover, we will see that within 
the Loudon-Fleury theory, there is a very strong enhancement of Aig intensity at small 
transferred frequencies cui — cuf = 5uj'^'^^ ~ 4( Ji 72)"'^/^, where Ui and cuf are the frequencies 
of the incoming and the scattered photon, respectively, and Ji is the intralayer exchange 
coupling. At the frequency shift Su^^^, the intensity in the Aig channel actually turns out 
to be larger than the intensity in the Big channel, which, as we recall, is nonzero already 
in the absence of J2. In the resonant regime, when the incident photon frequency becomes 
comparable to the single particle excitation gap of the insulator, there exists substantial A^ 
scattering already for a single layer. Nevertheless, we will argue that even in this situation, in 
a two-layer system, there is a measurable change of the Raman intensity near the frequency 
shift Su^^'^. We will also consider the scattering in Big geometry and will show that in this 
geometry the effect of J2 is much weaker than for Aig scattering. 

In principle, the presence of the new features in the Aig Raman scattering allows one 
to find the value of the interlayer coupling. In reality, however, in YBCO, the frequency 
shift 6uj^'^'^ is in the region where the dominant contribution to the Raman intensity comes 
from phonon rather than two-magnon scattering. However, we will argue that there are 
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still several features of the Raman profile in YBCO which are absent in the single-layer 
compound Sr2Cu02Cl2 and which allow to find an estimate for J2. We found from our 
analysis that J2 is likely to be about O.lJi. More rigorously, we can place the upper boundary 
for J2 as J2 ~ 0.25J. Neutron scattering data for the spin wave velocity Csw at half-filling 
| r?| yield Ji ~ 120meV. As Csw is only weakly dependent on J2 for small J2/ Ji, this implies 
that the probable value is J2 ~ 12meV, and the upper boundary for J2 is 30meV. This 
is consistent with the estimate bmeV < J2 < 20meV for J2 extracted from the analysis of 
NMR data on the double-layer material Y2BaiCu70i5 [|l9l, but substantially smaller 
than J2 ~ O.55J1 inferred from infrared transmission and reflection measurements and 
also substantially smaller than the theoretical estimate J2 = 56meV by Barriquand and 
Sawatzky ||21[ . 



The paper is organized as follows. We start in Sec. |I| with the Hubbard model at 
half-filling which has long range antiferromagnetic order in its ground state, and derive in 
a diagrammatic formalism the effective Loudon-Fleury model for Raman scattering in Aig 
and Big geometries in the non-resonant regime (i.e. assuming that the photon frequencies 
are smaller than the Mott-Hubbard gap). In Sec. |T|, we use this model to compute the 
Raman intensity in the Aig and Big channels first without magnon-magnon interaction. 



and then by including multiple scattering between magnons. In Sec. |rV| , we will discuss 
the resonant regime. Finally, in Sec. 0, we summarize our results and discuss them in the 
context of experimental data for single-layer Sr2Cu02Cl2 and double-layer YBa2CusOQ,i. 



II. THE DERIVATION OF THE LOUDON FLEURY HAMILTONIAN FOR THE 

NON RESONANT CASE 

In this section we derive the effective Loudon-Fleury model for Raman scattering using 
a momentum space diagrammatic formalism. This formalism has recently been applied 
to derive the Loudon-Fleury Hamiltonian for a single-layer system [0. Technically, the 
calculations for two-layer systems are more involved as one has to double the number of 
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fermionic operators. Conceptually, however, our approach is exactly the same as in Ref. ||12|| , 
and we therefore refrain from discussing the calculational steps in full length. 

The starting point of our calculations is the simplest one-band Hubbard Hamiltonian 
for a two-layer system at half-filling on a square lattice given by: 

H = -t J2 i'^hcj,a + dl^dj^^ + h.c.) -t'Y^ {cladi,a + h.c.) + UY^ ni^a,'\'ni,a,l (1) 

<i,j> i i,a 

where the c and d operators represent the electrons of layer 1 and 2, respectively, a = 1,2, 
^i,i,cr = cl^Cia, nj 2,cr = dl^di^a-, and t' is the hopping amplitude between the planes for which 
we assume t' < t (see Fig. |l]). We will also assume that t/U ^1 and thus perform our 
calculations only to leading order in t/U . This large U one-band model with only nearest- 
neighbor in-plane hopping is indeed a simplification, but it was argued in [|1^ that this 
model already contains the relevant physics for the analysis of Raman scattering in YBCO 
compounds. Here we follow this reasoning and assume that the model of Eq. (|1|) is valid. 

The mechanism of two-magnon Raman scattering is straightforward and has been dis- 
cussed a number of times in the literature P2|: the incoming photon with frequency uJi 



creates a virtual particle-hole pair, which then in turn emits two magnons with momenta 
q and — q before it annihilates into an outgoing photon with frequency uof. For the dia- 
grammatic calculation of the required matrix element for Raman scattering we thus need to 
compute two types of vertices: the vertices for the interaction between the electrons and the 
electromagnetic vector potential of the photons and the vertices for the interaction between 
the electrons and the magnons. 

The procedure to derive the coupling of light to the electrons was previously described by 



Shastry and Shraiman ||23|: the photons introduce a slowly varying vector potential A(r,t) 
in the presence of which the hopping term in the kinetic energy of the electrons acquires a 
phase {i-^ // A(I, t) -rfl^ . The Hubbard Hamiltonian is then expanded to leading orders in A. 
One further introduces the staggered magnetization as the spin density wave (SDW) order 
parameter and transforms to new fermionic operators which diagonalize the Hartree-Fock 



factorized Hamiltonian pj] . To derive the interaction vertex between fermions and magnons. 



one has to compute the transverse spin susceptibihties with and without momentum transfer 
Q = (vr,7r), and construct the interaction Hamiltonian which reproduces all dynamic spin 
susceptibihties. This procedure is unique though rather involved for a two-layer system. 
We will skip the details and present only our final results (some of the useful formulae are 
collected in Appendix |^). 

In terms of the new SDW quasiparticle operators, the Hubbard Hamiltonian takes the 
form 



(2) 



where 



and we set the lattice constant oq = 1. Here and below, the prime in the summation sign 
indicates that the summation is restricted to the magnetic Brillouin zone, i.e. to momenta 
where z/k > 0. For a two-layer system we obtain two pairs of conduction (described by the 
a and e) and valence {b and /) operators. The energy dispersions for each pair are shifted 
by Q = (7r,7r). The direct gap 2 A between the bands is determined by a self-consistency 
equation for the staggered magnetization and reduces to 2A = U in the strong coupling 
limit. 

The interaction between the fermionic current and the vector potential of light, which 
is relevant for two-magnon scattering at transferred photon frequencies small compared to 
the Mott-Hubbard gap, has the form 

^^.■ = -^Ej,'A-,. (4) 

jq is the current operator whose momentum can be safely set to zero since the velocity of 
light is several orders of magnitude larger than the Fermi velocity. To lowest order m t/U 
the components of the current operator are given by 

9ek 



k.cr " 



+ eJca/k.a + k.C. 



(5) 



We see that to lowest order m t/U the interaction with hght only leads to excitations of 
quasiparticles between the valence and conduction bands of each pair. Excitations within 
each band are of higher order mt/U and are therefore neglected in our present strong cou- 
pling theory. It is noteworthy that the fermionic current rewritten in terms of quasiparticles 
which decouple the Hamiltonian into two separate terms (i.e., a and (3 operators introduced 
in (|A3|) ) , contains only densities of these fermionic quasiparticles. As a result, there are no 
terms in jq=o which would correspond to excitations between the valence band of one pair 
and the conduction band of the other pair (e.g. aj^ af^o) HI • 



We now present the result for the Hamiltonian which describes the magnon-fermion 
interaction. A systematic way to derive this Hamiltonian is to extend the Hubbard model 



to a large number of n = 25 orbitals at a given site, and use a l/S" expansion p6|j27[| . The 
resulting spin-wave spectrum is that for a spin S antiferromagnet. In this Section, we will 
consider only noninteracting spin waves, i.e. we will keep only the leading linear term in the 
1/ S expansion. To simplify notations we will, however, not keep the overall factors of S in 
the formulae and thus present the results for n = 2S = 1. 

The spin-wave excitation spectrum of the two-layer antiferromagnet consists of one 
doubly degenerate branch with the dispersion 



fii(q) = 2Ji^(l - i/q)(l + i/q + J2/2J1) (6) 

for momenta q in the first Brillouin zone. The dynamic transverse spin susceptibility has 
poles aX Q = ^^i(q) and Q = fii(q + Q). For a single-layer antiferromagnet, ^^i(q) = 
f^i(q+ Q), and the two poles are indistinguishable. For two-layer systems, however, f2i(q) 
and fii(q + Q) are different, and it is convenient to introduce two types of magnon operators 
mq and riq with the dispersions l^i(q) and ^2(^1) = ^i(ct + Q)- 

The electron-magnon interaction Hamiltonian can be obtained in the same way as for a 
single-layer antiferromagnet |T2[: it is uniquely defined by the requirement that it should 
reproduce the forms of the dynamic spin susceptibilities, both with zero momentum transfer 
and with momentum transfer Q. The susceptibilities are presented in Appendix 0, and the 
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interaction Hamiltonian which reproduces them has the form 



H, 



el— mag 



ak+q,<7«k -a"^lqr_+(k, q) + ej,_^q ^^Cfc _amLqr__ (k, q) 

+ ^k+q,aV-^"^-qr++(k + q) + /k+q,a/k,-<T?^-qr + _(k, q) + 4+q,aek,-<T^-q^-+ (k, q) 
+ 4+q,a«k,-an^_q^'— (k, q) + ,,6k,-a^^-q^ + - (k, q) + foj^+q ,,/k,-a^^-q ^'++ (k, q) 

+ («k+q,aV-^ + 4+q,a/k,-a)"^tq$-(k, q) + 6j,+q,^ak,-a + /il+q,^ek,_^)mLq$+ (k, q) 

+ (4+q,a/k,-a + ej,+q ,,6k,-<T)^^-q5_(k, q) + (/ll+q,^ak,-a + 4+q,aek,-a)^^-q5+ (k, q) 
+ /l.C. 



(7) 



The vertex functions in Eq. |^ are 



$±(k, q) = Ak ± d , H±(k, q) = v^Af^q ± e 



r-..-(k,q)^^ 



(k,q) 
(k,q) 
(k,q) 



1 

1 



1 



^q(ek+q + 4 ) + '7q(ek+q ' 4) 



^q(ek+q + ek ) + ^q(ek+q " ^k > 



■^ql^k+q + ^k ) + 'Cq(fk+q ^k . 



■Cql^k+q + ^k ) + ^q(^k+q ~ ^k , 



(8) 



where, e.g. the upper signs are for r__|_ and the lower ones for F , and 




^q 



_ / l + ;/q + J2/2J1 ' 
V2\ 1-^q , 

1 /l-Z/q+ J2/2J/ 



1/4 



1/4 



eq=^ 



1 + 



(9) 



As for the case of a single-layer, the vertices which involve fermions from both, conduction 
and valence bands are of order U whereas the vertices involving only valence or only con- 
duction band fermions are of order t and in the non-resonant regime can be omitted in the 
calculations to lowest order in t/U . 

We now have all necessary tools to calculate the Raman matrix element Mr in a dia- 
grammatic technique. A simple experimentation shows that, just as in the case of a single 
layer, there are three diagrams which contribute to Mr to leading order in t/U (see Fig. 



One has to keep in mind, though, that the calculation of these diagrams now involves four 
bands and two different type of magnons. The integration over the internal frequencies in 
the diagrams of Fig. ^ is straightforward, and we obtain: 



1. for the diagram in Fig. ^ 

Mr =±4 (Al,(q)+/.i,(q))lE'^^"^ 



k,a,fj \ ^ 



4A 



dk 



4A2 - 9? 



(10) 



where the upper sign is for M^^^^^, and the lower sign for M^^^, the indices 1 and 
2 indicate whether the final state contains two magnons of type rriq or of type riq, 
respectively, and is a frequency equal to uji or uif which are indistinguishable to 
leading order in t/f/; 



2. for the diagram in Fig. 



M 



(1,2) 
R2 



-4 Ai,2(q) /ii,2(q)^E'( 



dkf 



-ef/3 



8A 



8A(4A2 + Q^] 
4Airr^ + ~(4A2^^T]2jI" 



3. for the diagram in Fig. ^ 

=4Ai,2(q) /^i,2(q)^E'f 



dkn 



dkf 



8A(4A2 + fi2 

(4A2 - 02^2 



;i2) 



As before, ej j are the polarization vectors of the incident and scattered photons, respec- 
tively, and the coefficients 2(q) and Ai^2(q) are defined as 



y2/ii(q) = r/q + r/q, V2Ai(q) = r/q - r/q, 

V2/^2(q) = + ^q, V^A2(q) = - ^q ■ 



(13) 



Simple algebra yields the relations 

1/2 



/^i,2iqj 



1/ 4J1 + J2 

.2V2^^i,2(q) , 



Ai,2(q) 



4JiZ/q ± J2 



1 / 4Ji + J2 

2U^i,2(q) 



1/2 



(14) 



I 4JiZ/q ± J2 

Comparing these results with those in Ref. [|T^ we observe that the form of Eqs. (|T^)-(|T2p 
is the same for a single- and a double-layer system. The information about the coupling 
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between the layers is therefore only contained in the coherence factors /ii,2(q) and Ai,2(q)- 
This is a direct consequence of the fact that the photons' vector potential only couples to 
the in-plane fermionic current. 
Finally, we use the relations 



1 ^ '/9ek \ /<9ek-q 



AT 5Z 



.)(■ 



-T'i 



de-L 



^ix^fx COS ~\~ ^iy^fy COS Qy 
^ix^fx ~l~ ^iy^fy 



(15) 



substitute them into Eqs. 



M 



(1,2) 
R 



-8t' 



2A 



4A2 - ^2 



'|12D, and obtain for the total Raman matrix element 



2Ai,2(q)/ii,2(q)(eixe/x + eiyCfy) =F 
(-^1,2 (q) + Aii,2(q))(eix'e/x cosg^. + CiyCfy cos qy) 



(1,2) 



(16) 



Notice that Mi^\q) = M}f^(q + Q). Since for the calculation of the Raman intensity 
we have to integrate over the whole Brillouin zone we can restrict our consideration to only 
one type of magnons, and just multiply Mr by \/2. 

We now change tracks and compute the matrix element Mr within the Loudon-Fleury 
theory [|TB|, i.e. we assume that the spins interact via the Heisenberg Hamiltonian 



<i,j>,a i 



(17) 



As before, a = 1,2, and the scattering of hght is described by the Loudon-Fleury Hamilto- 



nian 



Hl-F — A ^ y^ix^fx'^j,a ■ ^j+&x,a + ^iy^fy^j,a ' ^j+Sy,, 



(18) 



Here, A is a coupling constant and 6 = {Sx,6y) is a vector to nearest neighbors sites in a 
plane. Observe that this scattering Hamiltonian has the same form as for a single layer. 
This is again a consequence of the fact that the light only couples to the in-plane fermionic 
current. Following the standard procedure, the spin operators are now transformed to 
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boson operators via the conventional Holstein-Primakoff or Dyson-Maleev expansions. The 
easiest way to proceed is to introduce just one Bose field with momentum in the full first 
Brillouin zone, and then perform a unitary rotation to magnon operators which diagonalize 
the quadratic form for a two-layer Heisenberg antiferromagnet. As for a single plane, the 
transformation to magnon operators involves the same coefficients Ai(q) and /ii(q) as in Eq. 

(0)- 

Retaining only the term in the scattering Hamiltonian which contains two magnon cre- 
ation operators, we obtain for the Loudon-Fleury matrix element 



Rlf 



-A[2Ai(q)yUi(q)(ei^e/a,. + eiyCfy) - 

Ai(q) + /U?(q))(eia;e/^ cosg^ + ei^e/y cosg^) 



(19) 



Comparing the two expressions for Mji, Eq. and Eq. |T9|, we see that they are identical 
provided we identify the coupling constant 

2A 



4A2 - fi2 



(20) 



This concludes our derivation of the Loudon-Fleury Hamiltonian for a two-layer system. 

Before we proceed with the calculations of the Raman intensity, we would like to comment 
on the form of the Loudon-Fleury Hamiltonian. In some phenomenological theories for a 
single layer, the interaction Hamiltonian between light and spin degrees of freedom is written 

as 



(21) 



where 



1 
2 



(22) 



This formula is obtained from (0) if the term which describes scattering in Aig geometry is 
neglected. For a single layer, this procedure is legitimate as the scattering Hamiltonian in 
Aig geometry commutes with the Heisenberg Hamiltonian, and consequently there is no Aig 
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scattering. However, for a two-layer system, more care is needed as the Heisenberg Hamilto- 
nian now contains an extra term with an interlayer couphng, which does not commute with 
Hlf- As a resuh, if we want to rewrite the Loudon-Fleury Hamiltonian for two layers using 
the projection operator P, we necessarily have to introduce an extra term which contains 
spins from two different planes. Specifically, the Loudon-Fleury Hamiltonian expressed in 
terms of projection operators should have the form 

Hlf = Abi E ^(^^^ ^j, 5)S^j ■ S^j+s + Aai E ^(e*, e/, 0)Sij ■ Saj (23) 

j 

where P{ei,ef,6) is the same as before. Here, the first term describes scattering in Big 
geometry, while the second term (which couples spins from different layers) contributes to 
Aig scattering. Comparing the two forms of the Loudon-Fleury Hamiltonians, we obtain 
for the two coupling constants 



2A 



4A2 - fi2 



■Ji 



2A 



4A2 - 9? 



(24) 



The two forms of the Loudon-Fleury Hamiltonian indeed yield the same Raman matrix 
element as in Eq. (p^Q]). 



III. THE RAMAN INTENSITY IN THE NON RESONANT REGIME 

A. The Non-Interacting Case 
The Raman scattering cross section is proportional to the Golden Rule transition rate 



^ = JS"" E I ^« I' ^i^^i - ^^f + - e/) (25) 

where i and / are the initial and final states of the system, e^j are the corresponding 
energies, and the summation runs over all possible initial and final electronic states. Let 
us first neglect final state magnon-magnon interaction. Then, — e/ = 2fii(q), and using 
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(0) we obtain for the Raman intensity in the Aig and Big channels (dropping an identical 
overall prefactor) 



lB,,m oc 2^ (mJj'))' 5(?iu:, - hiUf - 2ni(q 



= 2A^ (l + A)^^ j^^j s{hu:^. - hu, - 2nM)) (26) 

where Q = uji — ujf and we have defined = {cosqx — cosqy)/2. We see that the Raman 
intensity in the Aig channel is proportional to {J2/ JiY thus vanishes with vanishing 
J2. In Big geometry, the changes in 1^^^ imposed by the interlayer coupling are minor: the 
form of the matrix element is preserved and the only changes appear in the prefactor and in 
the magnon energy dispersion Vti (q) . The momentum sums in Eq. (pB]) can be conveniently 
reduced to complete elliptic integrals. The resulting expressions for the intensities are, 
however, rather involved; they are collected in Appendix y. Here we discuss only the main 
features of the Raman spectrum. 

Our key observation is the following: The magnon energy f2i(q) is gapless at the zone 
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center q = 0, and has a gap, i^i(Q) = 2^ J2J1 at q = Q. Then, for f2 = cjj — u;/ < 2f2i(Q), 
only magnons with momentum near q can be excited. The numerators in both scattering 
geometries vanish at q = 0, and it is not difficult to show that the contributions from the 
q ^ region yield oc (see also below). For Q > 2r2i(Q), however, also magnons 
with q Q can be excited. In Big geometry, the numerator in Ib^^ contains the factor iy^ 
which vanishes at q = Q such that the opening of a new scattering channel does not cause 
substantial changes in the intensity which still scales as fi^. However, in Aig geometry the 
numerator in I^^g at q = Q is just a positive constant. In this situation, the scattering 
intensity changes drastically at Q = 2Qi{Q): using Eq. (p6D, we find 



_ /A,,(2^]l(Q) + M - /A,,(2r]i(Q) - M _ 255 j,^ 

IaJ2Qi{Q) - 600) {^i{Q)y ■ ^ ^ 
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For J2/J1 = 0.1 we obtain AJ^e/ ~ 1600, i.e. the enhancement of the Aig signal at the 
threshold frequency is very strong. 

Moreover, the value of the intensity in the Aig channel right above the threshold, Ia^^^ oc 
J|/fii(Q) oc Ji{J2/Jif^\ has the same order of magnitude as the intensity in the Big 
channel, Ibi^ oc (i7i(Q))'^ oc {J2/ Ji)^^'^- We found that for all reasonable values for J2/J1 
the ratio of intensities is ~ 1.7— 1.9. In other words, if the Loudon-Fleury approximation is 
applicable, and if one claims to observe the two-magnon profile in Big geometry at around 
f^i(Q), one should also observe, in a two-layer system, the signal of an even larger intensity 
in the Aig geometry. 

The intensities for the Aig and Big geometries without final state interaction are plot- 
ted in Fig. ^ for two different values of J2/ Ji- There are unphysical singularities in both 
intensities at the maximum magnon energy, but just as in the case of a single layer, they are 
artifacts of neglecting interactions between magnons. We will see in Sec. [1IIB| that once an 
interaction is included, the unphysical singularities are removed. 

B. The Interacting Case 

We now analyze how the two-magnon profile changes when the interaction between 
magnons is included. Fist of all, the magnon-magnon interaction renormalizes the spin- 
wave spectrum. To leading order in 1/5, which we only consider here, this renormalization 
can be absorbed into the renormalization of the exchange couplings 

where the momentum sums run over the whole first Brillouin zone (notice that we defined 
Vti without a factor of 25*). This renormalization comes from one-loop diagrams (Oguchi 
corrections Beyond the leading order in l/S", one has to solve Eq. (^) self-consistently 



and also include corrections with higher number of loops. Numerically, however, it turns out 
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that the dominant correction, at least to order l/S*^, still comes from one-loop diagrams ^1 . 
In other words, the actual magnon dispersion nearly preserves the same form as in linear 
spin wave theory, but contains renormalized coupling constants Ji^2- Below we will assume 
that this renormalization is already included into the definitions of Ji,2 and neglect it in our 
further consideration. 

Strictly speaking, to justify this approximation for the calculation of the effects due to a 
final state magnon-magnon interaction, we also have to prove that the dominant renormal- 
ization of the four-magnon interaction vertex can be absorbed into the same renormalization 
of the exchange integrals. We did not perform explicit 1/5* calculations for the vertex. How- 
ever, as the vertex itself has a factor 1/S in comparison to the magnon frequency, and all 
calculations involving magnon-magnon scattering will be performed only to the leading or- 
der in 1/S, then whether to use a bare or a renormalized Ji_2 in the vertex is beyond the 
accuracy of our calculations. For simplicity, we will henceforth use renormalized values of 
Ji 2 everywhere. 

We now consider in detail the renormalization of the two-magnon profile due to multiple 
scattering of two magnons. The magnon-magnon vertices can be immediately obtained 
from the Heisenberg Hamiltonian by applying e.g. the Holstein-Primakoff transformation 
to boson operators and a subsequent canonical transformation to magnon operators which 
diagonalize the quadratic part of the spin wave Hamiltonian. A detailed study of the effects 
due to magnon-magnon interaction in a single-layer antiferromagnet was already performed 



by Canali and Girvin |jTT[, and we follow here their line of reasoning. To leading order in 1/S, 
we can restrict ourselves to the scattering process which conserves the number of magnons. 
The effective scattering Hamiltonian then takes the form 



mag— mag 



1 

N 



(29) 



k,q 



with 



V{k,q)=A 



^]l(k)^]l(q) 



q-k 



A' 



+ 1 



(30) 
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and 



A=(J2 + 4Ji)/2, 5k= (J2 + 4JiZ/k)/2 . (31) 



In terms of A and i?k the magnon dispersion is given by ili(k) = y — B^^. 

In order to find the full vertex function for repeated two-magnon scattering we need 
to sum an infinite series of ladder diagrams. In Big geometry, the "side" vertices from the 
electron-photon coupling scale as 7q = (cos qx — cos qy)/2, and it is easy to see that the only 
term in Eq. (|30|) which contributes to scattering is the one with z/q_k. The evaluation of the 
ladder diagram series then proceeds exactly in the same way as for a single-layer system 
| rT| , |32[| . The analytical solution is presented in Appendix 0. 

The plots for the Raman intensity Is^g for two values of J2/J1 are shown in Fig. ^. The 
unphysical singularity that we found in the non-interacting case disappears, as expected, 
and we observe a pronounced two-magnon peak. We see that with increasing J2/J1 the 
two-magnon peak not only shifts to higher frequencies but that the amplitude of the signal 
also slightly increases. The latter, however, is mainly due to the overall factor (1 + J2/4J1) in 
the matrix element and a renormalization of the Loudon-Fleury constant A from magnon- 
magnon interactions. The shift of the peak position towards higher frequencies can be 
understood in the simple picture that the incoming photon flips two neighboring spins on 
the same layer. This creates misaligned spin pairs and thereby increases the total energy 
of the system. Evaluating the corresponding energy increase for a Neel state, we obtain a 
two-magnon peak at Q = 3Ji(l + J2/3J1) which is roughly consistent with what we find. 

In Aig geometry, the solution of the ladder series is more difficult since the "side" vertex 
behaves as ~ (1 — ^/q)/^7l(q) where q is the magnon momentum. At small q, this vertex 
scales linearly with q as a consequence of the Adler principle ||3^: the Raman matrix element 



includes the interaction between fermions and Goldstone bosons, and this interaction should 
vanish at the points where the magnon energy turns to zero. Because of the extra power of 
momentum in M^, the Raman intensity without final state interaction scales as /^^^ oc 
at very low frequencies, as mentioned before. However, the form of the "side" vertex in Aig 
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geometry, is not reproduced at the magnon-magnon vertex, and we in fact have to solve 
a set of coupled integral equations in order to get the result for the full J^^^. The explicit 
expression for Ia^^ is rather cumbersome, so we present it in Appendix |^ and here discuss 
only the key features of the solution. 

As the "side" vertices for Aig are invariant under transformations of the symmetry group 
Di^h of the square lattice, we can restrict ourselves to only that part of the scattering potential 
V^(k, q), which has the same symmetry, i.e. 



l^(k,q) = -2Ji 



(32) 



4Ji 4JiV 4Ji/ fiq l^k 

We see that V(k,q) actually tends to a finite value for k = q = 0. The cubic frequency 
dependence of lAig {^) is therefore actually an artifact of neglecting the final state interaction. 
When this interaction is included, Jai^ scales linearly with Q at the lowest frequencies. We 
also found that the real part of the polarizability has a logarithmic singularity at = 
2^1 (Q). This singularity gives rise to two effects: first, it makes Iai^ a continuous function 
of frequency, in other words, eliminates a jump in the intensity at 2fii(Q). Second, it gives 
rise to a strong peak in Ia^^ at frequencies somewhat smaller than 2Qi{Q,). Specifically, we 
found that near 2f2i(Q), the dominant contribution to the intensity comes from the third 
term in Eq. (32), and J^i, has the form 



where 

remains finite at 1] = 2^7i(Q) and is proportional to J|/f2i(Q). We see that there 
exists a very narrow peak in Ia^^ located at Q = 2i7i(Q)(l — exp(— 27r^7i(Q)/J2))- The 
intensity right at the peak is very high, Ia^^ oc l/r2f(Q). At small J2, the peak position is 
exponentially close to 2Qi{Q,). However, at larger J2, we found numerically that the peak 
is actually located at frequencies significantly smaller than 2i7i(Q). This last result agrees 
with the calculations of the two-magnon absorption profile in Ref. [EO]. 



The solutions for lA^g for two different values of J2/ Ji are graphically presented in Fig. 
The Aig Raman spectra have been evaluated on a finite lattice with 1000 x 1000 lattice 
points. A finite imaginary part i5 has been added to the energy denominators of the spin 
wave propagators in Eqs. (p6|), (P^), and ( |D8| ) in Appendix 0. This allows to study the 
influence of damping on the Aig two-magnon spectra. Without damping, the imaginary 
part of the polarizability has a jump at = 2fii(Q), and by the Kramers-Kronig relation, 
the real part of polarizabiliy {Ri term in Eq. (|3^)) necessarily has a logarithmic singularity, 
see Fig. |^c. With damping, the singular behavior near the threshold frequency is removed 
and the peak position is shifted closer to 2r2i(Q). Note also that, as in Big geometry, the 
divergence at twice the maximum spin wave frequency is removed due to the final state 
magnon-magnon interaction. 

Finally, for the purpose of comparison with experiments, it is useful to compute the 
ratio of the Raman intensities for Aig and Big geometries right at their peak positions. We 
found that this ratio is actually very small: for J2I J\ = 0.1 it is about 0.009 whereas for 
J2/J1 = 0.3, it is 0.044. In other words, though the Aig intensity at the peak is larger than 
the intensity of the Big signal at the same frequency, the overall scale of the peak is only a 
few percent of the two-magnon peak in Big geometry. We therefore have to conclude that 
in the non-resonant regime where the Loudon-Fleury theory is applicable, the extra peak 
in Aig geometry can hardly be separated from the background signal. We now consider 
what happens in the resonant regime, i.e. when the incident photon frequency becomes 
comparable to the Mott-Hubbard gap. 

IV. THE RAMAN INTENSITY IN THE RESONANT REGIME 

In our derivation of the Loudon-Fleury Hamiltonian, we have chosen the diagrams to 
leading order in t/U under the assumption that the energy of the incoming and outgoing 
photons is much smaller than the energy gap between the conduction and valence bands. 
Under these conditions, all denominators in the diagrams Figs, ^-c were of order U which 
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in turn allowed us to omit all diagrams with intraband scattering. The situation, however, 
becomes different in the resonant regime where the energy of the incoming photon comes 
close to the Mott-Hubbard gap 2A ~ U. Actually, most of the experiments on two-magnon 
Raman scattering in parent high-Tc compounds have been done with visible light frequencies 
which are only 0{J) apart from 2 A. 



It was shown in Ref. [|12| that in the resonant regime, the diagrams with intraband 
scattering are more relevant than those which contribute to the Loudon-Fleury theory, and, 
moreover, the dominant contribution to Raman scattering comes from just one diagram 
shown schematically in Fig. This diagram yields a Raman matrix element = + 
M^, where 

/dek \/9ek-q \r ± \ ±1^ 

jVf± = 8i — y " ^ (35) 

^ki:/3 {uJi - 2E± + t5) {uj, - fii(q) - E± - E^_^ + t5) [ujf - 2E±_q + t5) ' 

One of the key consequences of considering the resonance regime is that there exists a 
nonzero signal in Aig geometry even for a single-layer. Indeed, the absence of the Aig signal 
in the Loudon-Fleury theory was related to a particular form of the interaction Hamiltonian 
Hlf which contained only spin degrees of freedom. The inclusion of the intraband processes 
modifies the form of the interaction Hamiltonian with light, in which case it no longer 
commutes with the Heisenberg Hamiltonian even when = ej = (x + y)/\/2. 

The feature of the diagram in Fig. |^ which makes it dominant for Big scattering in the 
resonance regime is that it allows all three denominators to vanish simultaneously leading 



to a triple resonance enhancement |T2|. For two-layer systems, we should check whether 
or not the rapid variation of I Aig near 2Qi{Q,), can be enhanced when the incident photon 
frequency is tuned right to the triple resonance value. We performed computations analogous 



to those in [0] and found that there is in fact no enhancement at the frequency threshold 
for Aig scattering because the occurrence of the triple resonance requires that the fermionic 
velocities at momenta ko [uJi = 2E^{ko)) and ko + q, where q is the magnon momentum, be 
antiparallel to each other. For q = Q, we evidently have Vk-Ek lko= Vk-^k Iko+Q; i-e. the 
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two velocities are parallel. In this situation, the integration over the fermionic momenta near 
ko gives zero because all poles lie in the same half-plane. We also performed more detailed 
calculations by expanding the denominators up to second order around ko. This actually 
makes the integral over k — ko finite, but still there is no singularity in at Ui — 2i?ko; so 
we do not expect any substantial enhancement of the Raman intensity in the Aig channel 
due to a triple resonance. 

Despite the absence of the enhancement, the diagram in Fig. ^ is still relevant in the res- 
onant regime simply because it contains three denominators which all are 0{J). Since there 
is no resonant enhancement, then, to first approximation, one can just set the denominator 



in (|35|) to a constant and consider the basic structure of M/j as imposed by the interaction 



vertices between magnons and fermions. Performing simple calculations, we obtained from 



At small frequencies, the contributions to the Raman intensity come only from magnon 
momenta near q = 0. We see from (|36| ) that for these momenta, scales linearly with the 



magnon momentum, just as we found in the Loudon-Fleury approximation. Clearly then, 
the full Raman intensity in the absence of magnon-magnon interaction scales as lAig oc 
at small frequencies. We studied the effects of the magnon-magnon interaction and found 
that, as before, the bare form of the side vertex is not reproduced in a perturbation theory 
for magnon-magnon scattering, and the finite state interaction gives rise to a linear, rather 
than cubic frequency dependence of Iai^- Moreover, as M^j^ does not contain J2 as the 
overall factor, it obviously gives a dominant contribution to lAig- This in turn implies that 
at Q < 2i7i(Q), the Raman intensity in a double-layer system is roughly half of the intensity 
in a one-layer system. 

For Q > 2r2i(Q), magnon momenta near q = Q also contribute to It is not 

difficult to verify that this extra contribution has the same dependence on J2/J1 as in the 
Loudon-Fleury theory. Accordingly, if we consider only M|j, we obtain qualitatively the 
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same form of the Aig intensity profile as in the Loudon-Fleury theory - the only difference 
is that now the intensities below and above the jump at 2i7i(Q) are of the same order of 
magnitude. Our result for the Aig intensity computed with Mjj' with final state interaction 
is presented in Fig. |^. Qualitatively, the intensity profile is the same as in the Loudon- 
Fleury approximation, but the new features are a substantial increase in the intensity above 
the threshold at 2fii(Q) and the flattening of the Aig intensity slightly above the threshold 
frequency. There may also be a very narrow peak slightly below the threshold frequency 
(just as we obtained in the Loudon-Fleury theory), which we do not see because of a limited 
numerical accuracy. In any event, however, the singular behavior at this peak is eliminated 
by damping. 

The total matrix element for the Aig scattering is a sum of Loudon-Fleury and triple 
resonance contributions. Without studying in detail the frequency dependence of the denom- 
inator in (|35|) we cannot compare the overall strength of and M^. In general, in the 
absence of the enhancement due to an actual triple resonance, the two contributions should 
have the same order of magnitude |T2|. Experimentally, however, the overall intensity (and, 
to some extent, the form) of the Aig Raman profile demonstrates a substantial dependence 
on the incident photon frequency. Besides, as we noted above, the Loudon-Fleury result for 
the Aig intensity at the threshold is more than 1000 times smaller than the Big intensity at 
its maximum, while the experimental intensity ratio is about 40 times smaller in the vicinity 
of the triple resonance in Big geometry, and even far smaller away from the resonance. It 
is therefore very likely that the Loudon-Fleury contribution to the Aig intensity is just a 
minor correction to the intensity given by the triple resonance diagram. Notice also that 
M^^ has exactly the same form as the second term in M^, and its inclusion will just change 
the relative strength of the two terms in (|36|). 
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V. DISCUSSION 



We first summarize our results. We considered in this paper two-magnon Raman scat- 
tering in a two-layer Hubbard model at half-filling. We applied the SDW formalism and 
derived diagrammatically the Loudon-Fleury Hamiltonian for the interaction between light 
and spin degrees of freedom. We found that in a two-layer system, the scattering in Aig 
geometry is finite already in the Loudon-Fleury approximation. Without final state inter- 
action, the intensity in this channel scales as fi^ at low frequencies. The magnon-magnon 
interaction effects are numerically small, but nevertheless they change the frequency depen- 
dence to lAig oc f2 at small frequencies. Furthermore, there is a very strong resonance near 
= 4(Ji J2)^/^, when a second scattering channel opens up. At resonance, the amplitude of 
the Aig signal is larger than the amplitude of the signal at the same frequency. We also 
argued that in the resonance regime relevant to experiments on parent high-Tc compounds, 
there is no enhancement of the peak intensity in Aig geometry due to the actual triple res- 
onance. Nevertheless, the diagram with three resonant denominators is dominant in this 
regime as it yields a finite Aig intensity even without J2. 

In Fig. ^ we present the experimental data for the Aig Raman intensity for the single- 
layer Sr2CuO-2Cl2 and the double-layer Y Ba2Cu-iOQ,i compounds |3^. We see that at 
transferred frequencies uo > 2000cm^^, the intensity profiles in the two compounds are 
similar. The sharp peaks at the low energy tail of the two-magnon band are due to resonant 
multi-phonon scattering that becomes strongly enhanced for excitations close to 2 A. Despite 
the overall similarity of the two figures, there are clear differences at low frequencies. The 
intensity in a single-layer compound continues to decrease at frequencies smaller than the 
resonance frequencies for phonon scattering, while the intensity for a two-layer compound 
flattens at frequencies somewhat larger than the resonance frequencies for phonon scattering, 
and remains flat down to the smallest measured frequencies. 

At the moment, we do not understand the origin of the background contribution to 
the scattering in YBCO, but it is unlikely that this background contribution is related to 
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scattering in a half-filled insulator. A more likely possibility is that the background is due 
to the fact that the measured YBCO compound has some finite amount of holes. In any 
event, however, we see that the intensity flattens at about 1800cm~^, and is rather flat at 
even lower frequencies. 

This behavior is consistent with our result for the resonant regime where the experi- 
ments have been performed: the Aig intensity evaluated for J2 = O.lJi flattens at about 
I.8J1 ~ 1800cm-\ and is roughly two times flatter at low frequencies than in a single-layer 
compound. Indeed, the theory also predicts that there should be a jump in the intensity at 
the threshold frequency. However, if we associate the onset of flattening with ISOOcm^^, we 
flnd that the jump occurs at about 1200cm~^, i.e., right at the frequencies where the Raman 
signal is presumably dominated by phonon scattering, so there are little chances to observe 
this jump directly. We therefore believe that J2 ~ O.lJi is a reasonable though indirect 
estimate of J2. 

Notice that J2 ~ O.lJi ~ ISOi^ is consistent with the estimate obtained from the analysis 
of the NMR data in similar systems. More rigorously, we can place the upper bound on 
possible values of J2 because whatever the interpretation of the low frequency measurements 
is, the data above 2000cm~^ clearly show no influence of the interlayer coupling. This in turn 
implies that in any event, the threshold frequency is lower than 2000cm^^, or J2 < 0.25Ji. 
Even this estimate is substantially smaller than J2 ~ 0.55Ji extracted from the data of 
infrared transmission and reflection measurements in YBCO [0. Given that inelastic 
neutron scattering measurements were unable to detect the optical spin wave branch in 
antiferromagnetic YBCO up to GOmeV, it was argued [0 that the gap, 2(Ji 72)^''^, should 
be larger than 60meV. We therefore conclude that J2 must be in the energy range 8meV < 
J2 < SOmeV. 
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APPENDIX A: DERIVATION OF THE ENERGY DISPERSIONS OF THE 
QUASIPARTICLE CONDUCTION AND VALENCE BANDS 

In this Appendix we derive the dispersions for the valence and conduction fermions in 
the one-band double-layer Hubbard model at half-filling. After Fourier transformation to 
momentum space the Hubbard Hamiltonian Eq. (|l|) takes the form 

k,cr ' ' k ' 

+ ^ E (4'+q,aCk,-<xC-k',-aCk+q,a + 4'+q,a4,-<xC?-k',-a4+q,a) • (Al) 
k,k',q,(7 

The presence of long range antiferromagnetic SDW order in the ground state implies that 



T-r E (4+.,i^k,i) = -Tf E (4+.,T^k,T) = m . (A2) 



k k 



Introducing the linear combinations 



ttk.a = "^("^ka + dka) , /3ka = y|(^ka - C^ka) (A3) 



and decoupling the interaction term with the expectation values of (|A2|) , the Hubbard 
Hamiltonian (lAlD turns into 



H = ^E'I (-4it'k - t')("k,a«k,<x - Pl+^^„Pk+n,a) " U m sgn((T) \al^Pk+^,a + Pl^^^^ak, 



k,(7 



+ (4ti/k - ^')("k+7r,<7«k+7r,<x - /5k,a/^k,<T) - U msgn(a) [ak+^,^/?k,a + /^k.^ak+Tr.a 

where the primed momentum sum is restricted to the magnetic Brillouin zone. Two 
separate Bogolyubov transformations applied to the first and the second part of the 
Hamiltonian (|A4|) yield two pairs of conduction and valence bands with the dispersions 
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±E^ = ± J(4tz/k + t'y + A2 and ±E^ = ±J{4tui, - f )^ + A^- The self-consistency condi- 



tion for A = Um requires that 



1 
U 



2N 4" 



1 1 

+ 



EC 



(A5) 



APPENDIX B: TRANSVERSE SUSCEPTIBILITIES IN A DOUBLE LAYER 

ANTIFERROMAGNET 



The dynamic, transverse spin susceptibihty is obtained from the time ordered correlation 
function 



(Bl) 



where the indices a, (3 = 1,2 denote the layer. In terms of fermion operators the spin raising 
and lowering operators S^^^ are expressed through 



(B2) 



where a are the Pauli matrices and the c and d operators describe the electrons from layer 
1 and 2, respectively. Summing the RPA ladder diagram series for the transverse suscep- 



tibilities as described in detail in the literature for a single-layer system P^J3^J35[| leads to 
the following results (for S=l/2) 

Jl(l-^^q) Jl(l-Z/q+ J2/2J1) 



Xii (q,q + Q,^) 



1 



UJ 



(B3) 



when the spins are from the same layer {xti = X22 ^ind 

Jl(l-^^q) , Jl(l-Z/q+ J2/2Jl^ 



xi2 iq,q>^j 



(q,q + Q,^) 



UJ 



UJ 



uj^-Ql{q)+i6 uj^-Ql{q)+i6 



(B4) 
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when the spins are from different layers {xt'f = X2i~)- The poles of the susceptibilities are 
at fii(q) and fi2(q) = ^i(q + Q), as they should. Using these results, one can construct 
the effective Hamiltonian for the magnon-fermion interaction Eq. (|^) by the requirement 
that it reproduces the forms of all these susceptibilities. 



APPENDIX C: RAMAN INTENSITIES WITHOUT MAGNON MAGNON 

INTERACTION 



In this Appendix we present the closed forms of lAig and Isig without final state magnon- 
magnon interaction. The actual calculations have been performed to leading order in 1/S. 
Here, we present the results for S = 1/2. We introduce the short notation Q = Q/2Ji and 



'16t 



2\2 r 



1 + 



2A 



J2 



1 + 



(16t 



2\2 



1 - U_ 



1 



2A 



n 2 



4A2 - ^]2 
1 



A 

4Ji 



2fi 



1 + J2/4Ji)2 - nyA 



where 



z/+,_ = -(J2/4J1) ± V(l + ^2/4Ji))2 - fiV4 . 
For Aig scattering geometry we then obtain: 



i) for < < 2 JJ2/J1 



iAAn) = agin) 



1 + 1^4 



-K{U 



where K is the complete elliptic integral of the first kind; 



ii) for 2JJ2/J1 <n< 2JI + J2/2J1 



iAAn) = agin) 



1 + iy+ 



iii) for 2JI + J2/2J1 <n<2{l + J2/U1) 



lA,m = a g{n) (1 - u+)K{t+) + (1 - v^)K{t^ 



(CI) 



(C2) 



(C3) 



(C4) 



For Big scattering geometry we obtain 



i) for < < 2 JJ2/J1 
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(C5) 



where E is the complete eUiptic integral of the second kind; 



ii) for 2y'J2/Ji < < 2y'l + J2/2J1 

Injn) = b g{n) 1(1 + \k{U) - Bit 



iii) for 2^1 + J2/2J1 < ^] < 2(1 + J2/4J1) 

/B,^(fi) = 6^7(^]) |(l-z/+)[K(t+)-E(t+) 



(C6) 



(C7) 



APPENDIX D: RAMAN INTENSITY WITH MAGNON MAGNON 

INTERACTION 



In this Appendix we outline our calculations of the full Raman intensity with final state 
interaction. Considering repeated two-magnon scattering we sum the corresponding series 
of ladder diagrams (see, e.g. [|ll|])- The resulting integral equation for the full vertex function 
reduces to a set of algebraic equations which allows for an explicit solution. We skip the 
details and list here only the results. 

For the Raman intensity in Aig geometry we obtain 



lAigi^) oc Im- J 



1 + jain) 



where 

7=^.(4j,)^(i+j|: 

Here we defined 

A{n) 
H(n) 

E{Q) 



(n) = p(2)(^) + 2 Ji [i^W(^])E(^]) + —p^^\n)A{n)] . (D2) 



a 



-p«(^]) + i^(l)(^])i7(^]) 



1 + iJ2[pW(^^) - i^w(fi)/^(^^) 

2JiPW(f]) 



1 + 2j,Q{n) 
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(Dl) 



(D3) 

(D4) 
(D5) 



and also 



du ( 1 — Ua 



G'o(q, VL + uj) G'o(-q, -uj) 



1 ^(l-v^ 



N 



q ,l]i(q) ) ^]-2^]l(q) + ^5 
/ 

q 



(D6) 



(D7) 



where Gq{(\,uj) = {uj — ^^i(q) + i5)~^ is the non-interacting spin wave propagator. 
For Big geometry we find 



Jb (^]) oc Im- 



1 + 4JiL(o)(fi) + I6JKJ2 + 4Ji)2[L(o)(fi)L(2)(fi) - + 



4Ji 



(D9) 



where 



L^'^\^) = ^E/^ 7^ Go(q,fi + ^) Go(-q,-^) 



27r (q) 



1^ 



^r(q) l^-2fii(q)+i5 



(DIO) 



This form is similar to the result for a single layer. The plots of the full intensities are 
presented in Figs. | and ||. 
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FIGURES 

FIG. 1. The system under consideration is a two layer antiferromagnet with intralayer exchange 
coupling Ji = 4t^/C/ and interlayer exchange coupling J2 = 4t'^/C/. 

FIG. 2. The diagrams which contribute to lowest order in t/U to the Raman matrix element 
Mr. The fermions from the valence (conduction) bands are denoted by a dashed (solid) line. The 
emitted magnons are represented by a solid wavy line whereas the incoming (wj) and outgoing (cjj) 
photons are given by dashed wavy lines. 

FIG. 3. The Raman intensities in Aig and Big geometry obtained in the Loudon-Fleury theory 
neglecting final state magnon-magnon interactions. The jump in the intensity in the Aig geometry 
occurs at the frequency 2i7i(Q) = 4(J2 Ji)^/^. The overall shape of the intensities is shifted towards 
higher frequencies with increasing J2/ Ji- 

FIG. 4. The Loudon-Fleury Raman intensity in Big geometry with final state interaction for 
two different values of J2I J\- 

FIG. 5. The Loudon-Fleury Raman intensity in Aig geometry with final state interaction. 
The figures (a), for J2I J\ = 0.2, and (b), for J2I J\ = 0.4, include effects of magnon damping which 
were modeled by adding a finite imaginary part ib to the energy denominators of the spin wave 
propagators. The transferred frequency in these two figures is given in units of the maximum spin 
wave frequency ^max- The low frequency peak in the intensity is located at a frequency somewhat 
smaller than 2J7i(Q). With increasing damping the peak frequency is gradually shifted closer to 
2r2i(Q). For comparison, in Fig. (c) we plotted the intensity for J2/J1 = 0.4 without any magnon 
damping. 

FIG. 6. The "triple-resonance" diagram which gives the dominant contribution to the Raman 
intensity in the resonant regime. The notations are the same as in Fig.^ 
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FIG. 7. The "triple-resonance" diagram contribution to the Raman intensity in Aig geometry. 
The final state interaction is included, and J2/J1 = 0.1. Observe the flattening of the intensity 
above the threshold. There should be (in the absence of damping) a real jump at the threshold 
frequency - its smearing in the figure is due to limited numerical accuracy. The dashed line 
represents the intensity in the absence of J2. 

FIG. 8. Experimental Raman scattering data in x'x' scattering geometry for (a) single-layer 
Sr2Cu02Cl2 and (b) double-layer YBa2Cu30Q,i single crystals. The data are taken from Ref. 



|36|. Observe the flattening of the intensity in YBCO at around 1800cm~^. For both crystals, 
the excitation energy uJi is 2A + 2.9Ji (the actual values are 2.33 and 2.09 eV correspondingly). 
The continuum intensity at high frequencies is presumably due to multi-magnon Raman scattering, 
the sharp peaks at low energies are due to resonant multi-phonon scattering that becomes strongly 
enhanced for excitations close to 2A. The dashed line is a fit to linear -|- cubic frequency dependence 
/ oc [c(w/Ji) + (w/Ji)^], where c = 1.6 for Sr2Cu02Cl2 and 1.3 for Y Ba2Cu'j,0fi,i. 
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